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Generalized Morse Wavelets 
as a Superfamily of Analytic Wavelets 

Jonathan M. Lilly and Sofia C. Olhede 



Abstract — The generalized Morse wavelets are shown to constitute 
a superfamily that essentially encompasses all other commonly used 
analytic wavelets, subsuming eight apparently distinct types of analysis 
filters into a single common form. This superfamily of analytic wavelets 
provides a framework for systematically investigating wavelet suitability 
for various appUcations. In addition to a parameter controlling the time- 
domain duration or Fourier-domain bandwidth, the wavelet shape with 
fixed bandwidth may be modified by varying a second parameter, called 7. 
For integer values of 7, the most symmetric, most nearly Gaussian, and 
generally most time-frequency concentrated member of the superfamily 
is found to occur for 7 = 3. These wavelets, known as "Airy wavelets," 
capture the essential idea of popular Morlet wavelet, while avoiding its 
deficiencies. They may be recommended as an ideal starting point for 
general purpose use. 



I. Introduction 

Analytic wavelets — complex-valued time/frequency localized fil- 
ters with vanishing support on negative frequencies — provide the 
basis for a powerful analysis of oscillatory signals, see e.g. |1| and 
references therein. In order to study time- and/or frequency-localized 
variability of a signal x{t), the continuous wavelet transform 



W^{t,s) 



■r 



x{t) dr 



(1) 



is formed. This is seen as a set of bandpass operations with the scale 
noiTnalization ri = 1, or as a set of projections with the normalization 
n = 1/2, indexed by the scale parameter s. 

A practical question is which wavelet to use and why. In addition to 
the popular but only approximately analytic Morlet wavelet, a variety 
of analytic wavelets have been proposed, including the Cauchy- 
Klauder-Morse-Paul, Derivative of Gaussian, lognormal or log Gabor, 
Shannon, and Bessel wavelets |2|-(11- However, this profusion of 
wavelet types serves as a barrier to acquiring the knowledge necessaiy 
to carry out practical applications. In the absence of a suitable 
unifying theoiy for wavelet behaviors, the choice of a particular 
wavelet for a particular problem may even appear arbitrary. As 
attempts to address this, one may find in the literature very many 
examples of papers — frequently with overlapping content — in which 
different types of wavelets are applied to the same signal, in order to 
decide which wavelet is "better" in some sense. The large effort spent 
on such case-by-case analyses underscores the need for a systematic 
and unified framework for understanding analytic wavelets and their 
properties. 

Here we show that the generalized Morse wavelets ||6)-||9) effec- 
tively unify all wavelet types mentioned above, as well as the analytic 
filter and in a sense the complex exponentials themselves, as special 
cases within a single very broad family. These wavelets are defined 
in the frequency domain as 



*I'/3,-y(a;) = / ipis^^(t)e ^'^^dt = U{oj) ap^^, u'^ e 



(2) 
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where a^,-, is a normalization constant, U{ui) is the unit step function, 
and /3 and 7 are two parameters controlling the wavelet form[j The 
decision of what wavelet to use is greatly simplified by examining the 
roles of /? and 7 in controlling wavelet properties, as we undertake 
here, building on recent work [.9 J. 

The generalized Morse wavelets were introduced by Daubechies 
and Paul f6l as the eigenfunctions of a time/frequency localization 
operator, later examined in detail by [8 1. Their name reflects the fact 
that with 7 = 1, these wavelets become equivalent to a solution to 
the Schrodinger equation examined by Morse 1.10 J . Their frequency- 
domain form is essentially that of the probability distribution known 
as a generalized gamma or Stacy distribution fll], [12]. This distribu- 
tion has been rediscovered and re-examined multiple times [13], and 
finds applications ranging from droplet physics 1131 to econometrics 
|T4| to communications 1151 , 1161 . Consequently, the new results 
herein and in (9) regarding the generalized Morse wavelets are also 
relevant to the use of these functions as probability distributions. 

' The generalized Morse wavelet have an additional degree of freedom, the 
order k, which permits the construction of families of increasingly oscillatory 
wavelets orthogonal to a chosen wavelet of the form (|2), see e.g. 18]. Such 
higher-order wavelets can be expressed as sums of wavelets the form J2) with 
different /3 and 7 values. This degree of freedom, which is primarily useful 
in multiple-transform methods [e.g. ,7J, will not be pursued here. 
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Fig. 1 . Parameter space for the generalized Morse wavelet superfamily, with 
the Heisenberg area shaded over /3 and 7 on log-log axes. Note that the 
Heisenberg area is undefined for /3 < 1/2. Thin solid lines are contours of 
the Heisenberg area from 0.5 to 0.55 with a spacing of 0.005. In this log-log 
space, curves of constant duration P^-^ = \/~0^ are diagonal lines. Values of 
Pfj.'-f of 1/3, 1, 3, 9, and 27 ai'e marked by the thick gray lines. The dotted 
vertical lines mark the 7 = 1, 7 = 2, and 7 = 3 wavelet families. The 
heavy solid line is the right-hand border to the localization region of |8], 
with the left-hand border being at 7 = 1. A thin black line divides wavelets 
with positive skewness in the frequency domain (on the left) from those with 
negative skewness (on the right). Circles mark the locations of wavelets shown 
in Fig. |2] Letters mark the locations of other wavelet or analytic filter families, 
with arrows denoting limits, as discussed in the text: "L" for the lognormal 
wavelets, "C" for the Cauchy wavelets, "G" for the Derivative of Gaussian 
wavelets, "A" for the Airy wavelets, "e" for complex exponentials, "S" for the 
Shannon wavelet, "a" for the analytic filter, and "B" for the Bessel wavelet. 
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Fig. 2. Examples of the generalized Morse wavelets for 25 different (/3, 7) pairs, with time-domain wavelets ipi3.-y{t) on the left and frequency-domain 
wavelets ^^^(lj) on the right. In both panels, the parameter 7 increases from left to right, while the parameter /3 increases from bottom to top. Both 
parameters take on the values 1/3, 1, 3, 9, and 27, con'esponding to integer powers of 3. The duration Pfi,.y = \/Wf is constant along diagonal lines, as 
indicated by alternating white and gray backgrounds, and increases from the lower left to the upper right in both panels. For the time-domain wavelets, a 
thin solid line is the real part of the wavelet, a thin dashed line is the imaginaiy part, a heavy solid line is the modulus, a horizontal dotted line marks the 
zero axis, and vertical dotted lines mark itP^ .y. The time axis has been rescaled so that the duration P/3,-, is the same width in each panel, and the wavelet 
magnitudes have been rescaled to obtain the same values at time t = 0. For the frequency-domain wavelets, a heavy solid line is the wavelet, a thin solid line 
is the Gaussian approximation, and a thin dash-dotted line is the quartic approximation obtained later from (|6) by neglecting the residual term in fp,y{uj). A 
vertical dotted line marks the location of the peak frequency i^p^'y, which always occurs at the same position on account of a rescaling of the frequency axes. 



All software associated with this paper is freely distributed to the 
community as a part of a MATLAB toolbox, as described in the 
Appendix. 

II. Generalized Morse Wavelet Parameter Space 

By varying /3 and 7, the generalized Morse wavelets can take on 
a wide variety of forms, in addition to the rescaling by s which 
dilates or compresses a given wavelet. A map of generalized Morse 
wavelet parameter space is presented in Fig. [T] The dots mark the 
locations of the illustrative wavelets shown in Fig. |2] in which 
both time and frequency axes have been rescaled for presentational 
clarity. The original derivation of generalized Morse wavelets, see 
(Sj-fSl, utilized the fact that they are eigenvalues of a time/frequency 
localization operator for 7 > 1 and /3 > (7 — l)/2, which is the 
region to the right of the line 7 = 1 and to the left of the heavy black 
curve. As the generalized Morse wavelets are essentially a probability 
distribution in the frequency domain, they can be associated with a 
coefficient of skewness, see eqn. (42) of |9|. The thin black line in 
Fig.|2]is the location of vanishing frequency-domain skewness, with 
positive skewness to the left of this curve and negative skewness to 
the right. 

The shaded quantity in Fig. [T] is the Heisenberg area |3|, a con- 
venient and conventional measure of time/frequency concentration. 
The Heisenberg area of a function i/)(f) is defined as A^ = at Uoj, 
where at and a^^ are the time-domain and frequency-domain standard 
deviations, respectively; see e.g. (21) and (22) of f9] for definitions. 
Exact expressions for the generalized Morse wavelet standard devia- 
tions are given in |9|. The Heisenberg area Ap^^ of the generalized 



Morse wavelets decreases with increasing (3 for fixed 7. For /3 < i, 
the generalized Morse wavelet temporal standard deviation at is 
unbounded, as is the Heisenberg area. In the vicinity of the zero- 
skewness curve, there is a minimum of the Heisenberg area that 
approaches its theoretical lower bound for the greatest possible 
degree of concentration at A/3,-y = i. With increasing /3 this curve 
approaches the line 7 = 3, corresponding to an important family of 
wavelets identified by fW\ as arising from an inhomogeneous Airy 
function and therefore called the Airy wavelets. 

The wavelets in Fig.[2]become more oscillatory in the time domain, 
and therefore more narrow in the frequency domain, both from left 
to right as 7 increases, and from bottom to top as /3 increases. 
The frequency-domain wavelets obtain a maximum value at the peak 
frequency uj/j^^ = {P/'yY ' , determined by the frequency at which 
the derivative of >I',3^^(a;) with respect to ui vanishes, see 19|. A key 
parameter is the rescaled second derivative of the frequency-domain 
wavelet evaluated at its peak frequency, 



,8,7 



Pi- 



(3) 



*/3,7(^/3,7) 

It is shown by |9| that Pfs^^/n is the number of oscillations at the 
peak frequency that fit within the central window of the time-domain 
wavelet, and therefore Pft,^ is called the wavelet duration. Pp^^ is 
constant along diagonal lines sloping from the upper left the lower 
right in Figs. [T] and |2] so the number of oscillations encompassed by 
the wavelets along these lines is essentially the same; for example, 
along the main diagonal in Fig. |2l Pp.-, = 3 and so the number of 
oscillations in the central window is Pp.-Jn ~ 1. 
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Variations of the wavelet form for fixed duration Pfj,^ — y/^ can 
be expressed in several ways. From (|2), we see that the parameter /? 
controls the behavior near zero frequency, while 7 controls the high- 
frequency decay. It can be shown \9} that the wavelet time decay is 
l''/'/3,7(^)/V'/3, 7(0)1 ~ 1/i , thus P and 7 control the time-domain 
and frequency-domain decay respectively. Therefore moving along 
curves of constant Pg,^ from the upper left to the lower right of Fig.|2] 
exchanges a strong long-time decay rate, from /3, for a strong high- 
frequency decay rate, from 7. Conversely, moving diagonally from 
the lower left to upper right, the wavelets becomes more oscillatory 
in the time domain, and more narrow in the frequency domain, in 
such a way that the ratio of time-domain decay to frequency-domain 
decay /?/7 remains constant. It can be seen that the small /3 and 
large 7 wavelets in the lower right-hand comer are problematic on 
account of undesirable time-domain sidelobes and extreme frequency- 
domain asymmetry; this is also region within which the generalized 
Morse wavelet form (|2j is no longer the solution to the localization 
operator problem of [6], [8|, indicating the practical implications of 
the operator framework. 

Locations of other important wavelets families and analysis func- 
tions are indicated in Fig. [T] by letters. Note that herein the constant 
ap^ry in ^ is chosen such that ^^^^(cj/j,^) — 2, as is sensible for the 
1/s scale normalization in ([T) that we prefer, see |l9l. It was shown by 
(9) that the generalized Morse wavelets encompass two other popular 
families of analytic wavelets: the Morse, also known as the Cauchy 
or Klauder or Paul, wavelets for the choice 7 = 1, and the analytic 
"Derivative of Gaussian" (DoG) wavelets for 7 = 2. Furthermore, 
(9) shows that Q with (^,7) = (0,0) gives Vo,o(i) = S{t) + ^ 
or »l/o,o(i^) ~ 2U{uj), which is recognized as twice the analytic 
filter, while the complex exponentials are approached in a sense 
of convergence of moments as /? approaches infinity for fixed 7. 
These tendencies are apparent in Fig. |2] the lower left-hand comer in 
Fig. [2^ resembles an analytic delta-function in the time domain, while 
the upper right-hand comer in Fig. |2]) resembles a delta-function in 
frequency. 

In Section Unl we will show that with 7 = 3, the generalized Morse 
wavelets closely approximate a Gaussian while remaining exactly 
analytic, thus implementing the essential idea of the Morlet wavelet 
without suffering its drawbacks. In Section [TVl we will compare with 
three other analytic wavelets, the lognormal or log Gabor ||4], (5), 
Shannon |3|, and Bessel |2| wavelets. The lognormal wavelets are 
found to be recovered in the limit of vanishing 7 and fixed -P/3,7, 
that is, in diagonal motion to the upper left comer, while the Shannon 
wavelet is recovered in the opposite limit of motion to the lower right 
comer. The Bessel wavelet is found to be very closely approximated 
by the generalized Morse wavelet at (/3,7) = (22, 1/10). Altogether 
this is a fairly exhaustive list of commonly used analytic wavelets 
and analytic analysis functions, all of which are encompassed by the 
generalized Morse wavelet form (|2). 

III. Concentration AND Gaussianity 
A desirable characteristic for a wavelet would be to combine 
the special properties of a Gaussian — which is in a sense the most 
time/frequency concentrated function — with the constraints necessary 
for an analytic wavelet. This is not achievable; some compromise 
is necessary, see e.g. (21 p. 36-37. It will be shown that for 
integer choices of 7, the 7 = 3 family best achieves a high 
degree of concentration and Gaussianity, while remaining, like all 
the generalized Morse wavelets, exactly analytic. The Heisenberg 
area achieves its theoretical minimum value of one-half for a function 
which is a Gaussian envelope multiplying a complex exponential, but 
such a function is not a wavelet because it is not zero mean. One could 
say that the idea behind the Morlet wavelet is to introduce a correction 



in order to recover a zero mean function, hopefully preserving the 
desirable properties of the Gaussian in the process. In practice, this 
works well for oscillatory wavelet settings, but not for time-localized 
wavelet settings. 

The Morlet wavelet is defined in the time domain and the frequency 
domain as L2J, ^9J, respectively 



*i.(cj) 



ai, e 



■i(. 



(4) 
(5) 



where a^, is a normalizing constant. Note that the oscillation fre- 
quency u is not the same as the peak frequency at which V&^oj) 
vanishes; the deviation between these two frequencies grows as f 
become small and vanishes as v becomes large, see Appendix A of 
[9J. For time-localized wavelets with small values of v, [9J shows 
that the departure of the Morlet wavelet from analyticity — that is, 
its frequency-domain support on the negative frequencies — leads to 
interference terms that significantly degrade the resulting transform. 
Here we will show that the correction term also generates deviations 
from the Gaussian foiTn as well as reductions in the time/frequency 
localization as measured by the Heisenberg area. 

The inverse Heisenberg areas of the generalized Morse wavelets 
families with integer values of 7 from one to six are shown in 
Fig. [3^, together with the Heisenberg area of the Morlet wavelet. 
For large values of the wavelet duration, all wavelet families shown 
here approach the limiting value of concentration at A^,^. = | or 
1/Afj^ry = 2. For wavelets having small values of the duration, such as 
one would use to analyze rapidly vaiying signals, the Morlet wavelet 
concentration degrades more rapidly as the duration decreases than 
does the concentration of the 7 = 2, 3, or 4 generalized Morse 
wavelets. For highly time-localized settings, the Morlet wavelet 
become significantly less concentrated than these 7-families. 

For high degrees of concentration corresponding to A^p > 0.55, 
the generalized Morse wavelet family for integer 7 having the highest 
degree of time/frequency concentration for fixed choice of duration 
P^^-y is the Airy wavelet family at 7 = 3. This appears different from 
Fig. [T] which shows a minimum in the Heisenberg area somewhat 
to the right of 7 = 3, with the minimum Heisenberg area occurring 
in the vicinity of 7 = 4 for small values of j3 around /? ~ 2. The 
difference between these two perspectives is that different 7 families 
should be compared for a fixed value of P/3,^ and not a fixed value of 
/?. Taking this into account as in Fig. [3^ shows that the 7 = 3 family 
is the most concentrated family apart from for very small values of 
the duration, less than about P/3_^/7r ~ i, below which the 7 = 2 
and then the 7 = 1 families become the most concentrated. The 
Heisenberg areas of all generalized Morse wavelet families becomes 
unbounded at /3 = i, corresponding to Pp^^ = ^7/2, and thus the 
smallest possible wavelet duration having a finite area decreases with 
decreasing 7, as indicated by the vertical lines. 

A convenient way to express the variation in forms of the general- 
ized Morse wavelets as a function of /3 and 7 is via a Fourier-domain 
expansion about a Gaussian centered at the peak frequency oj^^^. 



*/3,7(^/3,7) 

:exp<!-ip|,, (^-^-1) +fp,,{u)} (6) 
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Fig. 3. (a) The inverse Heisenberg area l/Afjry for the generalized Morse 
wavelets with integer 7 from 1 to 6 as a function of wavelet duration Pg^,a& 
well as for the Morlet wavelet. Dashed vertical lines marked location at which 
P = 2 for 22ch choice of 7. where the Heisenberg area becomes infinite. 
The gray horizontal line corresponds to A^.-^ = 0.55. (b) The squared inner 
product between the wavelet and the Gaussian approximation (6). 



where the deviation from Gaussianity is controlled by the function 



//3,7(^) = ^g('>' 



3)-P/3,7 






— [(7- 

24 L^' 



- 3)^ + 2] Pf 



0,1 



1 



+ A*^'^(c.). 



(7) 



This expansion has been accomplished by Taylor-expanding the 
wavelet logarithm, and making use of expressions for the derivatives 
of In '^ p^~i{lij) from Appendix D of J9|. Here Avpl !^{u]) is a residual 
that is defined implicitly as including the remaining error terms in 
jSJl not explicitly represented in Q. 

Observe that in the Gaussian expansion ^, Pis^-, appears as the 
inverse of the standard deviation of the Gaussian, so that l/Pp^^ 
is a measure of the frequency-domain bandwidth, while 7 controls 
the degree of deviation from Gaussianity for a fixed choice of 
Pfi,-/. This suggests that P^,^ and 7 are the natural coordinates for 
the generalized Morse wavelets, with Pp^^, setting the number of 



oscillations and 7 controlling the shape. When 7 = 3, the lowest- 
order deviation from Gaussianity — the cubic term — vanishes, while 
the quartic term obtains its smallest magnitude. The cubic teiTn 
controls the local asymmetry of the wavelet about the peak frequency, 
while the quartic term describes a local narrowing or broadening of 
the wavelet peak in comparison with the Gaussian form; both are 
controlled by 7. Therefore in the sense of the local shape of the 
Fourier-domain wavelet about the peak frequency up^-y, the 7 = 3 
wavelets are the most symmetric and most nearly Gaussian members 
of the generalized Morse wavelet family for a fixed Pp^^ and any 
integer or non-integer choice of 7. 

In Fig.|2j3, the Fourier-domain wavelets are compared with approx- 
imations from keeping only the quadratic (Gaussian) term in ^, or 
with keeping up to and including the quartic term. In general, for 
highly localized wavelets the quartic and even the quadratic approxi- 
mations provide excellent approximations up to large distances from 
the peak frequency. However, for wavelets with a large Heisenberg 
area and hence poor time/frequency concentration, such as in the 
lower right-hand comer of Fig. |2ji, the low-order Taylor expansion 
^ is only accurate in the immediate vicinity of the peak frequency. 

A more direct comparison between the wavelets and the Gaussian 
form is shown in Fig. |3j). This is the inner product of the wavelet 
with the Gaussian form, specifically the quadratic approximation from 
{6]l neglecting the fp,j{ij) term, after normalizing both functions to 
unit energy. For the Morlet wavelet, we use the expressions for the 
peak frequency and duration given in Appendix A of |9|. Not only 
is the Morlet wavelet less concentrated than the generalized Morse 
wavelets, it is also less similar to a Gaussian, despite the intention 
behind its formation. While all the wavelets become increasingly 
more Gaussian as the duration increases, the 7-family for integer 7 
which is always the most similar to a Gaussian for particular choice 
of P/3,T, is the Airy wavelet family at 7 = 3. 

TV. Special Forms of Generalized Morse Wavelets 

Timiting foiTns of the generalized Morse wavelets can be found if 
we rescale the frequency argument lj by the peak frequency 1^/3, -y, 
defining a new frequency-domain wavelet 



$/3,^ (oj) = V^/J,^ {uJi3,-yUj) 



(8) 



such that the rescaled wavelet always obtains a maximum at fre- 
quency a; = 1. This becomes, from (|2j, 

<^p,^ (uj) = 2U{uj) ^'^e^/'Ti^-'-'') (9) 

after choosing a^.-y such that *I'/3,7 (1^/3,7) ~ 2. Setting /? = P^/7 
with P a positive constant, this becomes 



^p2/^^^ (cj) = 2U{uj)e 



(jlnuj + l-uj"')P''/j'' 



(10) 



making use of the fact that to — e 
tends to zero is 



The limit of the exponent as 7 



lim 



^■0 7^ 



7 In oj + 1 



where we have expanded uj'' 
quently it follows that 



1 + lna;7 + -7 hi ui + 0(7 ) 

= _ip2j^2^ (11) 

in powers of 7lnaj. Conse- 



^7lnt 



lim $p2/^,^H = 2f/(a;)e-^^'"^'(") 



(12) 



emerges as the small-7 limit of the generalized Morse wavelets with 
/3 chosen such that Pfj^^ — ^fp^ remains fixed. The wavelets are 
lognormal in the frequency domain, with 1/Pa,.y playing the role 
of the standard deviation, and present a highly asymmetric shape 
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about the maximum at cj = 1. Wavelets of the form ( I12l l have 
been proposed, apparently independently, by Pl and f5|, and appear 
to be in widespread use under the names "log Gabor wavelets" or 
"lognormal quadrature wavelets." 

In the opposite limit of infinite 7 with /? = P'^/7, moving towards 
the lower right comer of Fig. [T] another wavelet emerges. The uj^ = 
LO^ '"' term in (|9) tends to unity, while the exponential contains the 
indeterminate form aj^/7^ which tends to ^cj^^^ by rHospital's 
rule. We then find that ([9} becomes 



lim >I> 



py-, 



2 
0, 



lim 2U{Lj)e~ 

7 ^00 

<u] <1 
cj < 0,cj > 1 



2rect 



(13) 



where "rect" is the unit rectangle function. Thus the generalized 
Morse wavelet with 7 approaching infinity and /3 approaching zero 
becomes a bandpass filter with constant amplitude for frequencies uj 
between zero and one, and vanishing elsewhere. This limit is identical 
with the definition of the Shannon wavelet *I's(a;) I3j p. 246], 
corresponding a sine function in the time domain 



V's(t) = TTsinc [ — ] e 



(14) 



The slow (1/i) time decay of this wavelet, however, is an undesirable 
property that limits its usefulness. 

Together the Shannon wavelet and the lognormal wavelets are 
important because they delimit the most extreme properties that 
may be obtained for any choice of P/3,7. Note, however, that these 
limits are quite different. As seen in Fig. [T] the Heisenberg area 
becomes unbounded in the Shannon limit, but asymptotes to a fixed 
value for each P^,^ with decreasing 7 in the lognormal limit. In 
the Shannon limit the generalized Morse wavelets approach a fixed 
form for any choice of P/3,7, whereas in the opposite limit of the 
lognormal wavelets, the wavelet form remains a one-parameter family 
controlled by P^,^. In particular, it is perhaps surprising that the 
generalized Morse wavelets can both closely approximate a Gaussian 
in the frequency domain, at 7 = 3 as discussed earlier, as well as 
the lognormal form ( I12t . 

This means that the only exactly analytic wavelet in common use 
not included in the generalized Morse wavelet family is the Bessel 
wavelet |2. p. 32], defined by, with our choice of normalization 



i'Bit) = 



2e' 



7r\/l ^ it 



Ki 2^r~^ 



(15) 
(16) 



where Ki(-) is the first modified Bessel function of the second kind. 
It is evident by comparing ( II6t and (|2j that the Bessel wavelet 
does not formally lie within the generalized Morse wavelet family. 
Nevertheless, a choice of (/?, 7) may be found for which the Bessel 
wavelet is very closely approximated. A measure of the similarity 
between the two wavelets is the magnitude-squared inner product 



2 
a/3,. 



i'Bit) 



i>*fi,i^t) 



yIZ<, l^s Wl^ du ^JjZ^ |V';3,7(")I^ d-u- 



■dt 



(17) 



A numerical search over a range of (/3, 7) values shows that q| .^ 
obtains a maximum value of 0.9995 for (/3, 7) = (22, 1/10). At this 
parameter choice, the two time-domain wavelets are visually virtually 
indistinguishable. 



V. Discussion 

This correspondence has shown that a number of commonly used 
continuous analytic wavelets and analysis filters, rather than being 
intrinsically distinct, may all be regarded as special cases of a single 
unified family, the generalized Morse wavelets. This means that the 
wavelet transform itself may be regarded as a continuous function of 
not one but two parameters in addition to time and scale. The duration 
or inverse bandwidth Pp^ — \/^ sets the number of oscillations in 
the wavelet, while a shape parameter 7 shifts the transform between 
qualitatively different families with the duration held fixed. 

The decision of what wavelet to use is greatly simplified by noting 
the special properties of the Airy wavelet family at 7 = 3, previously 
identified by [9J. These are found to be the most Gaussian, most 
time/frequency concentrated, and most symmetric member of the 
generalized Morse wavelet family for any integer 7. We show that 
the Airy wavelets reflect the spirit behind the Morlet wavelet more so 
than even the Morlet wavelet itself, and consequently form a natural 
choice for a default analytic wavelet for general purpose applications. 
Particularly if time-localized settings are to be attempted, the Morlet 
wavelet should be avoided and the Airy wavelet used in its place, 
reinforcing the conclusion of |9|. 

Of course, real-valued or non-analytic wavelets may be useful, 
particularly in the analysis of discontinuities ||3], but for analyzing 
oscillatory phenomena analytic wavelets are commonly favored. 
There may be occasions when other members of the generalized 
Morse wavelet family than the Airy wavelets would be more suitable. 
One instance arises when the signal of interest may have a form 
corresponding to one of the wavelet families; this is particular true of 
the Gaussian and Cauchy families which may provide useful models 
for time-localized or "impulsive" events that themselves resemble 
wavelets. Another common scenario is that oscillatory variability 
may exist close to the Nyquist frequency, in which case one would 
like to sacrifice symmetry of the wavelet for a more rapid high- 
frequency decay. Further investigation of the generalized Morse 
wavelet parameter space may yield useful information regarding 
the properties and possible uses of analytic wavelets, in particular, 
the relationship between the roles of the parameters /? and 7 as 
differentiation and "warping" operators explored by ||9j, and the 
variations of the wavelet foiTn. 

Appendix 

All software associated with this paper is distributed as a 
part of a freely available MATLAB toolbox called JLAB, avail- 
able at http://www.jmlilly.net The generalized Morse wavelets 
are implemented with morsewave, while various properties 
are computed in the functions morseprops, morsefreq, 
morsederiv, morsespace, and morsebox. The wavelet trans- 
form wavetrans uses the generalized Morse wavelets by default. 
Finally, makef igs_superf amily generates all figures in this 
paper. 

References 

[1] J. M. Lilly and S. C. Olhede, "On the analytic wavelet transform," IEEE 
Trans. Inf. Theoiy, vol. 56, no. 8, pp. 4135^156, 2010. 

[2] M. Holschneider, Wavelets: an analysis tool. Oxford: Oxford University 
Press, 1995. 

[3] S . yisAlaX, A wavelet tour of .signal proces.sing, 2nd edition. New York: 
Academic Press, 1999. 

[4] H. Knutsson, C.-F. Westin, and G. Granlund, "Local multiscale fre- 
quency and bandwidth estimation," in Proceedings of IEEE International 
Conference on Image Processing. Austin, Texas: IEEE, November 1994, 
pp. 36^0. 

[5] D. J. Field, "Relations between the statistics of natural images and the 
response properties of cortical cells," J. Opt. Soc. Am. A, vol. 4, no. 12, 
pp. 2379-2394, 1987. 



IEEE TRANSACTIONS ON SIGNAL PROCESSING 



[6] I. Daubechies and T. Paul, "Time-frequency localisation operators: a 

geometric phase space approach II. The use of dilations and translations." 

Inverse Prohl, vol. 4, pp. 661-80, 1988. 
[7] M. Bayram and R. Baraniuk, Nonlinear and Nonstationary Signal 

Processing. Cambridge University Press, 2000, ch. Multiple window 

time-varying spectrum estimation, pp. 292-316. 
[8] S. C. Olhede and A. T. Walden, "Generalized Morse wavelets," IEEE 

Trans. Signal Process., vol. 50, no. 11, pp. 2661-2670, 2002. 
[9] J. M. Lilly and S. C. Olhede, "Higher-order properties of analytic 

wavelets," IEEE Trans. Signal Process., vol. 57, no. 1, pp. 146-160, 

2009. 
[10] P. Morse, "Diatomic molecules according to the wave mechanics II. 

Vibrational levels," Physical Review, vol. 34, pp. 57-64, 1929. 
[11] E. W. Stacy, "A generalization of the gamma distribution," Ann. Math. 

Stat., vol. 33, no. 3, pp. 1187-1192, 1962. 
[12] J. H. Lienhard and P. L. Meyer, "A physical basis for the generalized 

gamma distribution," Quart. Appl. Math., vol. 25, no. 3, pp. 330-334, 

1967. 
[13] S. Hegyi, "A powerful generalization of the NED suggested by Peter 

Carruthers," in VIII International Workshop on Multiparticle Production. 

Singapore: Worid Scientific Publishing Co. Ltd., 1999, pp. 272-286. 
[14] M. J. Zhang, J. R. Russell, and R. S. Tsay, "A nonlinear autoregressive 

conditional duration model with applications to financial transaction 

data," J. Econometrics, vol. 104, no. 2001, pp. 179-207, 2001. 
[15] N. C. Sagias, G. K. Karagiannidis, P. T. Mathiopoulos, and T. A. 

Tsiftsis, "On the performance analysis of equal-gain diversity receivers 

over generalized gamma fading channels," IEEE Trans. Wirel. Commun., 

vol. 5, no. 10, pp. 2967-2975, 2006. 
[16] M. D. Yacoub, "The a-^ distribution: a physical fading model for the 

Stacy distribution," IEEE Trans. Veh. Technol, vol. 56, no. 1, pp. 27-34, 

2007. 



